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We derive a purely fermionic no-sea effective theory, featuring positive-energy states only for the 
Walecka model. In contrast to the so-called mean-field theory approach with the no-sea approxi- 
mation, where the Dirac sea is simply omitted from the outset, we turn to the relativistic Hartree 
approximation and explicitly construct a no-sea effective theory from the underlying quantum field 
theory. Several results obtained within these two approaches are confronted with each other. This 
sheds new light on the reliability of the mean-field theory with the no-sea approximation as well 
as the role of the Dirac sea. Restricting to 1+1 dimensions, we obtain new analytical insights into 
nonuniform nuclear matter. 
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I. INTRODUCTION 

The Walecka model, also known as (cr, uj) model or 
quantum hadrodynamics I (QHD-I), is a simple relativis- 
tic quantum field theory (QFT), mimicking the main fea- 
tures of the nucleon-nucleon interaction: Forces between 
nucleons (spinor field tp) ^re mediated by the exchange 
of mesons. As far as isospin symmetric matter (equal 
number of protons and neutrons) is concerned, the dom- 
inant one meson exchanges come from the scalar sigma 
(tr) meson and vector omega (w) meson. The Walecka 
model defined by the following Lagrangian [ll-Q 
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C4cr 



(1) 



where F'^" — d^tu'^ — d^uj^ is the kinetic term associ- 
ated with the vector meson field w^. The nuclcon field 
is assumed to have N flavors, i.e., '0V' = J2k=i etc. 
In a phenomenological context one mostly works with 
N ~ 2, where the two indices denote the proton and 
neutron fields, respectively. From the viewpoint of renor- 
malized perturbation theory, the cubic and quartic scalar 
self-interactions have to be included in the Lagrangian, to 
make the model rcnormalizable in 3+1 dimensions. Only 
then the respective counter terms ^ and ^ a'^ arise 
systematically, when rescaling the bare fields in Eq. ([T]) 
and splitting the Lagrangian into a physical and a coun- 
terterm part (see Sec. Ill Ap in the usual way (cf. e.g. 
Ref. i). 

As the physical coupling constants arc assumed to be 
large in this model, perturbation theory is inapplica- 
ble. However, mean-field theory (MFT) with the no- 
sea approximation, the relativistic Hartree approxima- 
tion (RHA), and the Hartree- Fock (HF) approximation 
are known to yield physically sensible results. The basic 
idea behind these approximations is to reduce the under- 
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lying interacting QFT to the problem of a single parti- 
cle moving within the effective potential generated by all 
other particles. Here we turn to the simplest approach 
of this type, the relativistic Hartree approximation. In 
this context, one of the main problems in any practical 
calculation is posed by the question of how to account for 
the infinite number of negative-energy states constituting 
the Dirac sea. Due to this difficulty, early mean-field type 
calculations in attempts to base nuclear physics on field 
theoretical models like the Walecka model, were typically 
done without the sea. Similar ansatzc were also consid- 
ered for the Gross-Neveu (GN) model family in the realm 
of strong interaction physics. 

Recently a new method to derive a no-sea effective the- 
ory (NSET) with positive-energy states only has been in- 
troduced and validated within the framework of the GN 
model family in 1+1 dimensions, where various exact an- 
alytical solutions are known. Moreover, new analytical 
insights were obtained [H. In the course of the NSET 
approach the Dirac sea is integrated out explicitly. Most 
strikingly this gives rise to new, Dirac sea-induced inter- 
actions (between positive-energy states) not present in 
the original Lagrangian. The coupling constants and in- 
teraction terms of the resulting effective theory are well 
defined and finite. The NSET is free of any UV diver- 
gences. As the derivation of the NSET involves only stan- 
dard Feynman diagrams, the approach is quite general 
and can be easily extended to different QFTs amenable 
to a mean-field type approximation. In the present work, 
we apply this method to the Walecka model in 1+1 and 
3+1 dimensions. While the Walecka model was originally 
formulated in 3+1 dimensions [l|, there are also studies 
in 1+1 dimensions, where calculational efforts are less 
and even some exact solutions are known Our main 
interest is in the model in 3+1 dimensions. However, as 
from the viewpoint of the NSET, the models in 1+1 and 
3+1 dimensions can be treated in a similar way, and on 
the other hand, fully analytical calculations are feasible 
in 1+1 dimensions, we find it worthwhile to consider also 
the model in 1+1 dimensions here. 

This article is organized as follows. In Sec. lUwe con- 
struct the NSET for the Walecka model. In the course 
of this, renormalization conditions are specified. As it 
is conventionally done for the Walecka model, we uti- 
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lize counterterm renormalization. In Sec. IIIII a test of 
the NSET is performed by restricting to uniform nu- 
clear matter. Section IIVI focuses on nonuniform nuclear 
matter in the model in 1+1 dimensions. Following the 
NSET approach, explicit analytical calculations can be 
performed here. Our main focus will be on the question 
of identifying an expansion parameter that allows us to 
perform consistent and systematic approximations of the 
full problem. We intend to provide some first insights in 
the abilities of the NSET approach in the treatment of 
spatially nonuniform systems. Section |V] discusses some 
general observations and results obtained for the Walecka 
model in the NSET approach. In particular we contrast 
our findings with those obtained within MFT, were the 
Dirac sea is simply ignored from the outset. Several dif- 
ferences are stressed. We end with a short summary and 
conclusions in Sec. I VII 



II. 
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CONSTRUCTION OF NO-SEA EFFECTIVE 
THEORY 

A. Prerequisites 



In the course of counterterm renormalization, the orig- 
inal Lagrangian ([T]) is written as £ = £' + Cct, with 
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and the counterterm Lagrangian 
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While the tilded parameters and the fields in Eq. ([T]) are 
bare and unphysical, in Eq. ([2]) the analogous (untildcd) 
parameters are considered as physical and measurable. In 
this sense, the counterterm Lagrangian is not added. The 
original Lagrangian is rather split into two parts. C is the 
physical part; Cct contains the infinite but unobservable 
shifts between the bare and the physical parameters. As 
their choice is scale dependent, the physical parameters 
have to be specified by renormalization conditions, stat- 
ing under which circumstances they have to be measured. 
Subsequently, the counterterms are tuned such that the- 
ory and measurement agree under these conditions. Note 
that the individual counterterms can be unambiguously 
assigned to the (divergent) primitive diagrams depicted 
in Fig. [TJ Choosing the counterterms so as to cancel the 
divergences, the theory is rendered finite. In result, a di- 
vergent primitive Feynman diagrams never occurs on its 
own, but always in its renormalized, finite form. 

In order to remove the divergences in Figs.[T];,d in 3+1 
dimensions, the scalar interactions ^ and ^ cr^ in the 
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FIG. 1: Divergent primitive Feynman diagrams. In 3+1 di- 
mensions, diagrams a)-e) diverge. In 1+1 dimensions, dia- 
grams a), b) and e) diverge, while c) and d) are finite. Dashed 
line: free fermion propagator; dash-dotted line: scalar propa- 
gator; curly line: vector propagator. 



original Lagrangian ([T]) are mandatory. Without incor- 
porating them, the associated counterterms proportional 
to and would not be generated, when rescaling 
the fields in Eq. ([1} to their renormalized counterparts 
in Eq. ([2)), and the respective divergences could not be 
cured. We will however assume (cf. Ref. 0) that the 
physical couplings, cubic and quartic in the scalar field, 
vanish in the vacuum at zero external momentum trans- 
fer, i.e., C3 = C4 = 0. This choice actually corresponds to 
two renormalization conditions, fixing the explicit values 
of and 04 (see below). 

In 1+1 dimensions, the diagrams depicted in Figs.[T};,d 
are finite. Hence, the counterterm contributions ^ 
and ^ cr^ are not needed. Working in 1+1 dimensions, 
we therefore set C3 = £4 = from the beginning. Conse- 
quently, the respective counterterm contributions are not 
generated and as = a4 = 0. 

In the following, the divergent fermion loop diagrams 
are treated in dimensional regularization. Our conven- 
tions are as follows 
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where d — D + 1 — e, with e in £> + 1 space-time 
dimensions. 

While we choose and determine the values of the vari- 
ous counterterms as in the original work of Chin , we 
consider it as helpful to shortly repeat the basic ideas 
of their determination in this context. In particular as 
our reasoning is somewhat different here and tailored to 
the formulation of the NSET. Instead of assembling the 
counterterms in a separate counterterm Lagrangian (cf. 
Eq. ([3])) and just resorting to them, when calculating ob- 
servables, the counterterms are an integral ingredient of 
the NSET Lagrangian, thereby rendering the NSET free 
of any UV divergences. 

We first determine the nucleon self-energy in the vac- 
uum. The vacuum is assumed to be translationally in- 
variant. In general the full fermion self-energy "Eh (cf. 
Fig. [2]) has scalar (S) as well as vector (VJj) contributions 
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FIG. 2: Dyson equation for fermion propagator in RHA. Ob- 
viously the last two terms in the second line do not arise 
in 1+1 dimensions. Dashed line: free propagator; solid line: 
dressed propagator; shaded circle: self-energy; dash-dotted 
line: scalar propagator; curly line: vector propagator. 



However the vector contribution can be shown to van- 
ish in the vacuum, resulting in a purely scalar contribu- 
tion ^h\vs.c ~ ^'^ ~ Here we introduced the physical 
fermion mass in the vacuum M = S'jvac- A common 
choice 0, H is 



= ^ M = m, 



(6) 



i.e., the mass m occurring in the Lagrangian C ^ C + 
CcT already corresponds to the physical mass in the vac- 
uum. In this sense, all possible loop corrections are al- 
ready included and Yih vanishes in the vacuum. It is 
obvious from Fig. [2l that the renormalization condition 
Eq. ^ is fulfilled by fixing 5m such that the single prim- 
itive Feynman diagram Fig.[T]a,) vanishes at zero external 
momentum transfer 
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The values for the remaining counterterms are deter- 
mined in the next step. Therefore we turn to the con- 
struction of dressed meson field propagators in the vac- 
uum. In the presence of the cubic and quartic couplings 
of the scalar field this is rather complicated, but signif- 
icantly simpler without them. Let us for the moment 
turn to the Walecka model in 1+1 dimensions where 
C3 — ~ and also — ~ 0. For this choice, 
the Lagrangian is quadratic in the scalar as well as the 
vector field. Hence, both fields can be integrated out 
exactly, yielding a purely fermionic theory. We directly 
turn to the renormalized version, i.e., integrate out the 
renormalized meson fields. This yields 



C = ijj (i0 — TO — 6m) "0 



2 ^ (l + a^)n + TO2 -a2 
1 



X I -S,f - ^ 1 ip-ff,ip . (8) 



From Eq. ([5]) one can easily read off the (momentum 



FIG. 3: Resummation of scalar vacuum polarization graphs 
coupled via a propagators into a dressed scalar field propaga- 
tor. 



space) propagators, incorporating the appropriate coun- 
terterms, for the a and uj'^ fields. They are given by 



iA(p) = 



(1 + aa)p^ - ml + q;2 + i?7 



(9) 



and 
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-Sf + 



ml 
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respectively. While the counterterms a2 and can 
be associated with iA(p), a^j belongs to \DJ^{p). The 
dressed versions of the propagators © and ([TU)) in the 
vacuum are now easily constructed. One simply has to 
account for all possible vacuum ( " — " ) fermion loop in- 
sertions in the respective meson propagator. To simplify 
the notation, the unrenormalized " — " loop with two 1 
insertions and external momentum transfer k is denoted 
by iE(A:2). We evaluate it by standard Feynman rules 
and expand the result in powers of the four-momentum 
k flowing through the graph. Let us first consider the 
scalar field propagator. In RHA it emerges as depicted 
in Fig. 131 Summing up the geometric series, the corre- 
sponding analytical expression is 



iA'(fc) = 



(1 



)fc2 



a2 - gl^ik"^ 



IT] 



(11) 



Note that both a2 as well as can be determined from 
A' (A:). In the following we turn to a particular choice of 
the counterterms. Interpreting the parameter to^ to be 
the renormalized, physical mass of the a meson in the 
vacuum, we demand that the effective scalar propagator 
is of the form 



iA'(A:) 



fc2- 



177 



+ 0(fc^ 



(12) 



In order to determine the effective vector field propaga- 
tor, one has to invoke the same steps as for the effective 
scalar propagator. We just have to replace the scalar cou- 
plings by vector couplings and the scalar propagators by 
vector propagators in Fig. |31 Due to the Lorcntz indices. 
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a) 



b) 



FIG. 4: Diagrams giving rise to the physical interactions a) ~ 
and b) ~ ct^. The divergent fermion loops are renormalized 
by an adequate choice of the counterterms as and 04- 



the determination of the effective vector propagator is 
somewhat more involved. Performing the resummation 
in terms of a geometric series, a matrix in Lorentz in- 
dices has to be inverted. To simphfy the notation, the 
unrenormahzed " — " loop with 7^ as well as 7^^ inser- 
tion at external momentum transfer k is referred to as 
in^'^(A;) in the following. Extracting the tensor structure 
in the Lorentz indices 



in^'^(fc) = iig'^'k^ - fc''fc'^)n(fc2) , 



(13) 



the vector propagator in the effective theory can be cast 
into the following algebraic form 
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X {-6/ + C{k')k,¥') + 



C{k^)k^ 
kM' 



(14) 



where we have introduced the momentum-dependent 
scalar function 



11] 



[55in(fc2) + ia, 



(15) 



Note that for C{k'^) = 0, the vector field propagator of 
the original theory, Eq. ([TU]), is recovered. We fix by 
implementing this choice. 

Let us now reconsider the situation. To simplify the 
construction of the dressed meson propagators we turned 
to 1+1 dimensions and omitted the scalar interaction 
terms ~ and ~ cr^ for a moment. However, it turns 
out that the expressions for the dressed meson propa- 
gators remain valid in 3-1-1 dimensions for zero physical 
couplings C3 = C4 = 0. As noted above, this choice cor- 
responds to two renormalization conditions fixing both 
and 0(4. The physical and cr* interactions are de- 
picted in Fig. m We fix 013 and in such a way that 
they exactly cancel the vacuum loops with three and four 
scalar insertions. Hence, with the choice C3 = C4 = any 
additional contributions to the dressed meson propaga- 
tors iA'(fc) and i£'^^(fc) in the vacuum are ruled out 
by the mutual cancellation of the diagrams depicted in 
Figs. [T^) and[T|D), respectively. Note that nonvanishing 
physical coefficients C3 and C4 would significantly com- 
plicate calculations, but are not required for a consistent 



model. Finally all counterterm parameters have been 
fixed. While 5m and the counterterms multiplying pow- 
ers of the meson fields (a2 • • • cti) are infinite, those multi- 
plying derivatives of the meson fields (ao-, ctu) are finite. 
The role of the latter is to ensure that the meson prop- 
agators have the expected momentum dependence. We 
are now prepared to construct the no-sea effective theory 
with positive-energy states only. 



B. No-sea effective Lagrangian 

Let us now turn to the problem of finite fermion num- 
ber. The RHA still has the same basic structure as in 
the vacuum, except that the nucleon self-energy is in gen- 
eral space-dependent and that the vector contribution in 
Fig. [2] no longer vanishes. The nucleon propagator gets 
an extra contribution from the positive-energy valence 
states. We denote the vacuum and valence particle con- 
tributions by "— " and respectively. 

Somewhat analogous to the decomposition of the nu- 
cleon Hartree propagator at finite Fermi momentum kf, 
featuring a mass M = M{kf) and accounting for a con- 
stant vector field V = V''{kf) 



iG«(p)=iGf(p)+iGg(p), 
iGf(p) 



\Gl{p) 



where 



E{p) 



5{f - E{p))6{kf - IpI) , (16) 



and E{p) ^ + A/^ we define 



iG^(p)=iG_(p)+iG+(p) 
iG_(p) 



|i — m + i?7 



(17) 



Note that in contrast to iG|? (p), which has the form of the 
fermion Hartree propagator in the vacuum, - but contains 
the mass determined at finite Fermi momentum M{kf), 
iG_ (p) corresponds directly to the fermion Hartree prop- 
agator in the vacuum [mass m = M{kf = 0))]. The 
second term iG^{p) can be associated with the positive- 
energy contribution. Its explicit form is not important 
for the following discussion. Hence, the one-loop dia- 
gram contributing to the self-energy for instance, then 
naturally splits up into the two pieces shown in Fig. [S] 
In the effective theory with the Dirac sea integrated out, 
diagram Fig. [5^) is accounted for by the mass term for 
the nucleon field 



(1) 



(18) 



(From now on the superscript on Cos refers to the power 
in the nucleon condensates ijjip and -ipj^ip of the respec- 
tive contribution to the NSET Lagrangian.) All the nu- 
cleon spinors appearing in the remainder of this article 
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FIG. 5: Decomposition of the one-loop contribution to the 
nucleon self-energy into a) negative and b) positive-energy 
contributions. 



refer to the sector, hence we will omit the "4-" label. 
Diagram Fig. corresporrds to the lowest-order no-sea 
RHA calculation with the interaction mediated by the 
original scalar propagator descending from Eq. ([1]). Ob- 
viously the dressed meson propagators can now be rein- 
terpreted as the effective propagators mediating between 
nucleon fields. They already contain all possible "— " 
loop insertions. Hence the four-fermion interactions of 
the effective Lagrangian can be written as 



C 



(2) 



(19) 



The presence of the propagators in Eq. ([T9| shows that we 
are dealing with nonlocal effective interactions between 
positive-energy fermions. Here, however, our aim is to de- 
rive a local no-sea effective theory, valid in the vicinity of 
the fully occupied Dirac sea only. We restrict ourselves to 
small external two-momentum transfer k in the effective 
propagators Eqs. ((TT|) and ([T?|) . This will enable us to 
derive an "almost local" effective Lagrangian well suited 
for analytical studies. Let us therefore trade the effective 
propagators for effective, fc-dcpendent couplings. As both 
the effective scalar field propagator and the effective vec- 
tor field propagator describe massive fields, the resulting 
couplings are well defined in the infrared and contain no 
singularities. We assume that the mass terms dominate 
the inverse effective meson propagators, i.e., rn^ 3> fc^ 
and 3> fc^, respectively. To keep the derivation of 
the NSET compact and general, we introduce coefficients 
Cijk in the following. As will become obvious below, the 
indices successively refer to the powers of (iptp), (V'T^V') 
and the derivative of the respective interaction terms 
in the NSET Lagrangian. Their explicit values in 1-1-1 
and 3-1-1 dimensions are listed in ((28| and (|29)) . While 
the effective scalar coupling is then given by 



mi ^ J 



(20) 



the effective vector coupling is determined by 



(5!,,cff(fc))/ 



mi 

X [l+C024fc'] +0(fc') 



(21) 



As a result, the nonlocal, effective interactions in Eq. 
decompose into an infinite number of local interactions 
(involving increasing powers of derivatives in position 
space, as fc'' — ;> —19^). We obtain 



C 



eff 



9a 



2ml 



9a 



2mi 



(Di^^p) (^V) 



2mi 



C204 



(22) 



and 



C 



(2,1-) _ 
off ~ 



9: 



2ml 

ji 

2m^ 

9l 
2m; 



('/'7i'V')(V'7''V') 

4 [{d.i^Yi^yd.i^Y^P) + {Ui,iui>){4>Yi')\ 



-C024 [(9.V^7»(n5^^7''V') 



(23) 



for the scalar and vector contributions, respectively. The 
ellipses denote higher derivative terms. Moreover, effec- 
tive interactions involving higher powers in the nucleon 
condensates are generated in the process of integrating 
out the Dirac sea. Note that, having derived the effec- 
tive couplings and (HH, the remaining steps in the 
construction of the no-sea effective theory for the Walecka 
model coincide with those invoked in the context of the 
GN model family Q. The only difference being that in- 
stead of scalar and pseudoscalar insertions, we now have 
to deal with scalar and vector, i.e., 1 and 7^ insertions in 
the "— " loops. In order to determine the effective interac- 
tion terms it is convenient to evaluate nucleon self-energy 
contributions and derive the effective interactions there- 
from. The respective self-energy contributions arc made 
up of internal "— " loops with more than two insertions 
and outermost loops, coupled via the effective point 
couplings defined in Eqs. ([20| and (|2T|) . In this article 
we will consider terms in the effective Lagrangian that 
are ~ (nucleon condensates)^ at most. Figure [6] depicts 
the self-energy diagrams giving rise to effective interac- 
tion terms of the power 3 [[Fig. [5^], 4 [Fig. [S})) and[Sl;)] 
and 5 [Fig. [6]i)-|6f )] in the nucleon condensates. To de- 
termine a given effective interaction term in the NSET 
Lagrangian, one has to sum up all topologically distinct 
diagrams with the respective number of outermost 
loops. Due to the infinite series making up the effective 
couplings, every of the diagrams in Fig. |6] represents an 
infinite series of loop diagrams. Remember that the mass 
term already includes an infinite number of loops due to 
the self-consistency condition (see Fig. [2]). Moreover it 
can be shown, that all "— " loops with at least one 7'^ in- 
sertion (and an arbitrary number of 1 insertions) vanish 
at zero external momentum transfer. In this article we 
do not explicitly derive all the interaction terms in the 
effective Lagrangian. We just illustrate the basic proce- 
dure for the effective interaction arising from Fig. [5^), 
evaluated with scalar interactions only. The underlying 
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FIG. 6: Contributions to the nucleon self-energy that are the 
seed for new interaction terms in the effective Lagrangian. For 
the sake of clarity, the "+" and "— " labels are only depicted 
in the first few diagrams and omitted in the remainder. 

idea is an expansion under the loop integral in powers 
of the momentum transfer. The analytical expression for 
the respective self-energy contribution is 

Sm= -i / ^^(27r)J(fc-fci-fc2) 

X (iff'.cff (fci + fc2)) {igl^ffih)) (i5',cff(fc2)) 

x{^ij)kAi^i^)k,, (24) 



where we treat both fcf and k!^ as of 0{k). Making an 
ansatz for the interaction term in the effective Lagrangian 



+ cMo^^)i^i^f + ■ ■ ■] (25) 



and deriving the corresponding self-energy contribution 
therefrom, the parameters C300 and C302 in Eq. (|25p can 
be inferred. The other effective interactions are deter- 
mined analogously. It turns out that there is just one 
further nonvanishing contribution from Fig. [S^, namely 
that arising from the "— " loop with 7'^, 7"^, and 1 inser- 
tion. In this case, the interaction term in the effective 
Lagrangian is given by 

- (^7,^)9" (^7/3 V)] +••• • (26) 

Note that due to the reason given above no new non- 
derivative interactions involving higher powers of the con- 
densate ■ip-f'^^ are induced. We do not continue in a fully 
systematic way beyond this level. However, in view of 
the applications we have in mind, we single out specific 
terms. Namely the non-derivative terms ^ (V'V')^ ^-nd 
^ {ipip)^. Collecting the different contributions, we fi- 
nally arrive at the following NSET Lagrangian 



222 2 

£efi' = - m)^p + ^{i>i^f - ^(□#)(V5V') + ^c204(n'^V')(V5V') - 

2 2 

' [{d,^Yi>){d^{^ri,) + (□^57.^) (^7"^ V')] - TrTCo24 [(a,^7''V')(n9^^7'^V') + {a^^j.^K^bj^^ 



' ^ ci22(^^) K i^lo.^) dp ii'l"^) - {^Ic^) d" (^7/3^)] + cM^^f 

Tfl, 



2\'^ /2\^ /2\'^ 



mi I \ mi. I \ mf 



C302(nV'V')(V'^)' + cMnr + C5oo(^V')' + • • • • (27) 



I 

While the general structure of the effective Lagrangian is the different interaction terms differ. We obtain 
the same in and 3-1-1 dimensions, the coefficients of 

^ Ngl 1 Ngl 



C204 — — T + -TT^ J , Co24 



rrvi 1207rm* ' m},, 307rm'* 
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C122 



C400 



C500 
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N 
207rm3 



C300 



N 
67rm 



N 



N9I 



5Ngl 
irrni 



C302 



27rm?. 



(28) 



in 1+1 dimensions and 



C204 



C122 



1 

N 



C024 



127r2m 



, C300 — , C302 



1 

N 



C400 — , C500 — - 



TV 



(29) 



in 3+1 dimensions. Let us emphasize again that the effec- 
tive Lagrangian of the Walecka in 3+1 dimensions does 
not include non-derivative interactions cubic and quar- 
tic in the fermion bihnear ipip. This is a consequence of 
our renormahzation conditions (cf. Fig.|4]). The determi- 
nation of the (V)^)^-interaction is therefore very simple 
in 3+1 dimensions. The only contributing diagram is 

Fig. EH). 



III. TESTING THE NO-SEA EFFECTIVE 
THEORY 



Wc propose to test the NSET Lagrangian for a low- 
density system of fermions with Fermi momentum kf, 
assuming unbroken translational invariance. Only in this 
(possibly unphysical) case do we have the exact analyt- 
ical solution to compare with. If the condensates are 
assumed to be spatially uniform, only the non-derivative 
terms enter and Eq. ([27]) reduces to the following effective 
Lagrangian 

4ff = ^(i^ - m)^b + ^i^^f - 

Zmt 2m, , 



m 



mi J 



(30) 



The Euler-Lagrange equation then allows us to identify 
the self-energy via 

^C',ff={i0-EH-m)i' = O. (31) 

Using the decomposition of the full self-energy into scalar 
and vector contributions, Eq. ([5|), Eq. ((3T|) becomes 



(i^ -M + 7^V;) ^ = . 



(32) 



The corresponding effective mass M is then given by 

M = m-^(^V)-3f^) C3oo(^V')' (33) 
mi \miJ 

By construction the condensates that appear in Eq. ([55)1 
only refer to the positive-energy sector. Using the famil- 
iar decomposition of the full fermion propagator (mass 
M) at finite Fermi momentum kj, into vacuum and va- 
lence particle contributions Eq. P7)) . the explicit expres- 
sion for the desired condensate can be easily inferred 



(2^)^ ^p2 + Af2 



(34) 



where D refers to the space dimension. This implies that 
the terms written explicitly in Eq. (|33p should allow the 
determination of the effective mass valid up to 0{kj) m 
1+1 dimensions and up to 0{kj^) in 3+1 dimensions, 
respectively. The lowest-order term omitted in Eq. ([55)1 
is ^ {ipil))^ . 

Note that the underlying assumption is that M can be 
expressed as a power series in kf. That this is possible is 
not a priori clear Q. However, here as well as in the case 
of the energy density, considered below, the exact RHA 
results arc known and confirm this assumption within the 
Walecka model. 

For illustration let us quote the explicit expression of 
the low-density expansion of M (expanding in fc^) in 3+1 
dimensions. It is given by 



M 



3 



10m2 



30 



144m6 
35 



1408to4 



3co 



16c2 



175m^ 



where 



Ng'^ 



(35) 



(36) 



Here, kf is assumed to be small as compared to the mass 
scales m and ma-. More precisely, the expansion is in the 
dimensionless ratios kf/m, kf/m^ and kf/m^, respec- 
tively. In contrast to the situation in the vacuum, the 
vector part of the self-energy also yields a finite contri- 
bution. As {ipfip) = for a spatially uniform system, it 
is determined as follows 



with 



PB = {tpH) = 



PbSq , 



(2^)^ 



(37) 



(38) 
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Let us now consider the exact solution. Demanding 
translational invariance, we determine the fermion mass 
at finite kf in 3+1 dimensions 



M - m 



Ngl 1 


[ Jo 


(27r)277i2 1 








+ In ( 


\ m j 



d^p 



M 



6mM'' 



3 



(39) 



The Taylor expansion of the exact result, Eq. (|39l) . repro- 
duces Eq. ([55]) . The expression for the vector self-energy 
associated with Eq. ([50]) . Eq. ([57)1 . agrees with the ex- 
pression for the full vector self-energy, determined under 
the assumption of translational invariance. This serves as 
a good test of our NSET in the special case of unbroken 
translational symmetry and confirms that the coefficients 
of the non-derivative terms in £off have been evaluated 
correctly in 3-1-1 dimensions. In 1-1-1 dimensions the cal- 
culation is performed analogously. 

No analytical bound state solutions are known for the 
Walecka model. Therefore we cannot test the derivative 
terms in the present case. The analysis of localized bound 
states in 3-1-1 dimensions would presumably require nu- 
merical methods and is outside the scope of the present 
work. However, we can obtain some analytical insights in 
1-t-l dimensions. In the following section we demonstrate 
that the treatment of nonuniform nuclear matter within 
the NSET approach is indeed possible. 



By a slight generalization of their work to the simulta- 
neous appearance of scalar and vector interactions, fea- 
turing two different coupling constants, this problem can 
be solved analytically. The derivation can be found in 
Appendix [X] This then enables us to organize and take 
into account higher order corrections to £cff in terms of 
a small parameter, at least in a certain regime. 

It turns out that there is a solution only for 



9o 



> 



9u 



(41) 



It remains to specify an expansion parameter. Therefore 
we invoke additional assumptions. Let us set 



1 



and 



1 



B 

2 ' 



(42) 



where A, B and 7 are dimensionless parameters and mo is 
a unit mass. Introducing three independent parameters 
A, B and 7, the problem is overdetermined. We can set 
one of the parameters A or B equal to one. This is done 
such that 



max 



1 1 



1 



7TO§ 



or, equivalently (cf. Eq. dH])), 

max (A, B) = 1 

and 

min (A, B) = C < 1 
Hence inequality (pij) translates into 



(43) 



(44) 



(45) 



IV. APPLICATION OF NO-SEA EFFECTIVE 
THEORY IN 1+1 DIMENSIONS 

So far, we have derived the no-sea effective theory for 
the Walecka model and tested it in a special case. Here 
we use our effective Lagrangian (|27p in 1+1 dimensions 
to study bound states of n < fermions. Our main fo- 
cus will be on the question of identifying an expansion 
parameter which allows us to perform consistent and sys- 
tematic approximations to the full problem. 

We start with the assumption that the leading-order 
interaction term in the effective theory is given by the 
two non-derivative four-fermion interactions 



— m)i/j - 



A. 

2ml 



2ml 



(40) 



In order to find a n-fermion bound-state solution associ- 
ated with the effective Lagrangian (|40)) . it is convenient 
to make use of an approach introduced by Lee et at. [loj 
to construct exact localized solutions of 1+1 dimensional 
field theories with four-fermion interactions. 



C < 



9l 



or C> 



(46) 



respectively. Moreover, recall the conditions ^ fc^ 
and ^ fc^, which were necessary in deriving a local 
effective theory. Independent of Eq. (|43)) . they translate 
into 7m^ > fc^ ^nd 7m^C"^ > fc^. As C < 1, both 
inequalities are simultaneously fulfilled if the single in- 
equality 7mg 3> fc^ is true. Here we want to use 7"^ 
as our (dimensionless) expansion parameter. This means 
that we also have to guarantee 7^1. The specifications 
on the expansion parameter can thus be summarized 

by 



77710 > and 7 > 1 



(47) 



Incorporating Eq. ()42p and the appropriate specifications 
for A and B in the exact ri-fermion bound-state solution 
of Eq. (PD|) , given in Appendix one finds 



1 



1 



~ — , ip^^ij) — , 



7 



(48) 



the last equation being a direct consequence of our choice 
of the 7 matrices. Consequently, we have 



1 



(^7»2 



1 



1 

7 



(49) 
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We are now in a position to count 7~^-powcrs of the vari- 
ous terms contributing to the fuh effective Lagrangian. In 
accordance with our assumption, interaction terms fea- 
turing higher powers of the fermion bihnears are of higher 
order in 7"^ and therefore subleading as compared to 
Eq. (|40| . Let us, for example, concentrate on the La- 
grangian vahd up to 0(7^^) 



9a 



N 

lA 

2mt 



9l 



2m^ 
3 



9Z 



2ml 
I9I 



2 mi 



(□^V)(V'V') 



(50) 



The fermion bound-state mass Af„ computed from this 
Lagrangian should be valid up to 0(7"^). As J dx ^ 7, 
the validity of the series expansion of M„ is reduced by 
one power of 7"^ as compared to the series expansion of 
£off- In order to determine M„. we follow the procedure 
developed above in the context of the NJL model and 
use a series expansion in 7^^. To do this, we turn to the 
equation of motion derived from Eq. (jSOp . written in the 
following form 



with 



S{x) = m + s ■ 



^ 9a 2 1 



27rTO 



(51) 



(52) 



and 



to 0(7 ^). For the self-consistent potentials, we obtain 



S^mil 



S22 



1 



S44 



1 



cosh^ z 7^ V cosh^ ^ cosh^ ^ 7 7^ 



-t- 



Vcosh^^ cosh"*^/ 7^ 
/ S62 , S64 , See 



(56) 



\ cosh^ ^ cosh ^ cosh^ C / 7 



and 



Vo 



[ P22 



1 



Pi2 



P44 



i cosh^ z 7^ 

P52 



1 

,cosh^^ ■ cosh^'^y 7^ 



+ 



cosh^^ cosh'^^y 7^ 
/ P62 , P64 , Pee 



(57) 



Vcosh^^ cosh"*^ cosh'^C/ T 



1 



The explicit expressions for the coefficients Smn , Pmn can 
be found in Appendix IB] In this particular example, we 
have demonstrated that a consistent truncation schema 
based on an expansion parameter can be found. As a 
result, we have obtained an analytically computable ex- 
pression for Mn in powers of 7"^. As pointed out above, 
exact bound-state solutions arc not known for this model 
so that we have derived new results. They could in prin- 
ciple be checked numerically. 

The analysis of localized bound states in 3-1-1 dimen- 
sions would presumably require numerical methods and 
is outside the scope of the present work. 

In the following section we rather take a step back 
and discuss several more general features of the NSET 
approach in the context of the Walecka model. In par- 
ticular we confront the NSET approach with the MET 
approach with the no-sea approximation. 



Vo (x) = po-\ 2" dl po 



(53) 



V. GENERAL OBSERVATIONS AND RESULTS 



expressed sclf-consistently through the condensates 



A. MET with the no-sea approximation 



S(x) = -4W> and po(x)=-4-(^^^)- (54) 
mi mt 



Here we used that (i/jj'^ip) vanishes [cf. Eq. (|48)) ]. We 
solve Eq. ([5T|) for a single energy level a = 0, assuming 
n < N. If Ag^ — Bg^ > 0, the bound-state mass is then 
given by 



^ 180 





1 


2A 




N 




TT 





^ ^Agl-Bgl 1 
1920 mg 74 



(55) 



2\ 3 



A2'Z T32'2 
o 2 ^ 9a- a 9u. 

Am c 2 



where wc defined c = [Ag"^ — Bg1)/mQ and kept terms 
up to order 1/7'^. A and B arc specified by Eq. and 
below. Note that the effective Lagrangian pO)) is valid up 



In order to clarify the differences, let us shortly re- 
view the MET approach in the context of the Walecka 
model 0|- MET can be seen as an ad hoc prescription 
to gain practical solutions for the complicated field the- 
oretic problem, posed by the Lagrangian ([1]). While the 
fermion field is still considered as a field operator, the me- 
son fields are replaced by their expectation values, which 
are assumed to be classical fields with finite values. The 
Dirac sea, occupied with an infinite number of negative- 
energy fermion states, is simply omitted. In result, there 
do not occur any divergences and the "bare parameters" 
in Eq. ^ are directly interpreted as their physical coun- 
terparts. As above, we do not allow for physical couplings 
cubic and quartic in the scalar meson field and hence as- 
sume C3 ~ ~ 0. In result, the Euler-Lagrangc equation 
for ip becomes 



[i7A.9^ - 9^1^.{^n - (m ^ga{a))]^ = 0. 



(58) 
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The classical meson fields satisfy the field equations 

(V2-m2)(a) =-g„(V;V), (59) 
{V^^ml){Lo^)=-g^{^^^). (60) 

Since one considers matter at rest, the classical three 
vector field (a;) is assumed to vanish Hence, the 
classical meson fields are generated self-consistently by 
the fermion condensates. 

We formally rewrite Eqs. ([5^ and (pO)) as 



1 



-9^ 



-{i>H), 



(61) 



and insert them in Eq. ([58|) . This yields the following 
equation of motion 



17^5'' + 7o5i 



1 



m + 5^ 



uj T-72 2 

1 



V2 -m 



= 0. 



(62) 



It is important to note that while the original field 
equations obviously refer to the full fermion field in 
MFT one uses them as equations for only, simply sub- 
stituting ij) for -0+. From a field theoretic point of view, 
such a drastic procedure seems hard to justify. However, 
the MFT approach turned out to describe many physical 
phenomena surprisingly well. 



B. Uniform nuclear matter 



Dirac sea is treated explicitly and enters the derivation 
of the effective theory in the form of " — " loops. This 
gives rise to inevitable divergences, whose treatment re- 
quires an adequate renormalization procedure. In the 
course of this, renormalization conditions have to be spec- 
ified. The resulting effective theory is dependent on these 
renormalization conditions. Their choice dictates the ex- 
plicit form of the various couplings in the NSET. We 
choose them such that the effective quantities resemble 
their bare counterparts in the vicinity of the fully occu- 
pied Dirac sea. Only then, having performed the lengthy 
derivation of the NSET, switching from bare to physi- 
cal parameters, the correspondence noted above can be 
established. In this respect the NSET approach indi- 
cates why (as long as a Hartrce treatment is justified) 
the MFT approach yields reasonable results for uniform 
nuclear matter at small Fermi momentum in the Walecka 
model. 

Note that the no-sea effective theory Lagrangian (PP]) 
should be valid up to C(fcj) in 1+1 dimensions and up 
to 0{kY) in 3-1-1 dimensions, the lowest-order correction 
term being ~ (i^fp)^ ~ kj and ^ {'4'4')^ 
lively (cf. the remarks in Sec. IIIip . 

It is easy to derive the energy of uniform matter at 
finite Fermi momentum within the framework of the 
NSET. One simply has to solve the single. particle equa- 
tions of motion, derived from the respective effective La- 
grangian for their energy eigenvalues and integrate 
them up. Note that in Hartree approximation this corre- 
sponds to integrating over the positive-energy solutions of 
the massive Dirac equation (mass Af), featuring a chem- 
ical potential Vq. In addition one has to account for dou- 
ble counting corrections (dec). One obtains 



kY , respec- 



Let us first focus on (uniform) nuclear matter. Keep- 
ing only the leading-order interaction terms in the NSET 
Lagrangian for spatially uniform matter, Eq. pop . we ob- 
tain Eq. (|40)) . The effective mass and the effective vector 
self-energy are then given by 



M = m 



V = — 



9o 



(63) 



and the associated equations of motion can be written as 



i2n) 



iEip)-Vo) +dcc. 



(65) 



where E{p) = ^/p'^ + ^P, with M from Eq. JM]). Spec- 
ifying to 3+1 dimensions, Eq. (j65p becomes 



2N 



' 6?p 
(2^ 

TV 



E{p) 



2ml 



2m2 



5m (27r)2 



to; 



(66) 



1(3— 771 



V' = . (64) 



Note that Eqs. ((63)) and (|64p arc the defining equations of 
the MFT approach to stationary, uniform nuclear matter, 
where (cr) and are assumed to be constants. Let us, 
however, emphasize that the MFT approach significantly 
differs from the NSET approach. 

In contrast to MFT, the NSET can be traced back 
to a relativistic Hartree approximation of the full, un- 
derlying relativistic quantum field theory, Eq. ([T|). The 



While the first three terms in Eq. (|66l) constitute the full 
MFT result, the last term written explicitly is the first 
nonvanishing NSET correction (~ fc}^)- Since the lowest- 
order term omitted in Eq. (j66p is proportional to {ijjijj)^ ~ 
k^ , the terms written explicitly should reproduce the 
exact result up to 0{ky). In this sense, the MFT result 
in 3+1 dimensions should be valid up to 0{kY)- 

These predictions can indeed be verified by turning to 
the series expansion of the exact result for the energy 
of uniform matter, obtained by performing the full RHA 
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calculation in 3+1 dimensions 



^RHA = 2N I ^ vV + M^ +77^(^1- 







(27r)3 



5^ 2 



TV 



2m 



2(27r)2 



2.9? 

, , M\ 25 
A/-" In — - — M 



m / 12 



+ AmM^ - 3m^M^ + -m^M - -m'^ !> . (67) 
3 41^' 



We now turn to the energy per nucleon, determined 
by £/pb- Obviously, the NSET result for £/pb with £ 
determined by Eq. (|66p is valid up to 0{kY) and can 
systematically be improved, while the MFT result can 
be trusted up to 0{kY) only. 



For illustration the series expansion of the NSET result 
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up to O(fci^) is stated explicitly 



PB 



Z{kff i{kff , {kff l^ikff , 2l{kfr 21(fc^) 
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539m7 lOOlm^ 



3 r 
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14000to5 



9lPB 



431200m7 

4 r 



rni AOtt'^ 10m 



mi. 8tt^ 



(68) 



Here, we take over the grouping, as well as the interpre- 
tation of the different contributions from Eq. (4.8). 
The first term is the fermion rest mass, followed by the 
nonrelativistic Fermi-gas energy and the first few rela- 
tivistic corrections. The next two terms (proportional to 
Pb) give the nonrelativistic limit of the potential energy 
coming from the vector and scalar mesons. The following 
term in brackets (with overall factor pb) is a rclativistic 
correction to the scalar potential energy that arises from 
the Lorentz contraction factor in the scalar density, eval- 
uated for fermions of mass to. The next three terms (with 
overall factors pg, p% and p%) are also corrections to the 
scalar potential energy. 

Let us shortly discuss the physical range of fc/, where 
such an expansion is sensible. Turning to the dimen- 
sionless quantity — , fc/ appears in the dimensionless 
ratios kf /m, kf /m^, and kf/m^ only. The smallest mass 
is governing the range of validity of the expansion. With 
TO = 939 MeV, m„ = 550 MeV, and m^ = 783 MeV 0, 
the smallest mass corresponds to m^ . Assuming the ratio 
to be small when it is < 0.5, we obtain the following esti- 
mate kf/m^r = fc//(2.8fm^^) < 0.5. Hence, we expect a 
good convergence for kj < 1.4 fm~^ Note that this is in 
agreement with the convergence behavior of the curves 
depicted in Fig. [7] (see below). 

The scries expansions of the full MFT and RHA re- 
sults for the energy per nucleon in uniform nuclear mat- 
ter agree with each other up to 0{kj^). By including 
higher-order interaction terms in the NSET Lagrangian, 
the NSET result can be systematically improved to re- 
produce the RHA result up to any desired order. The 
higher-order terms constituting the full MFT result can- 



not be justified in the NSET approach. 

To clarify the situation, we consider an illustration 
(cf. Fig. [7]) as helpful. However, several comments are 
in order here. Obviously, a plot of the energy per nu- 
cleon requires a choice of the coupling constants. Trying 
to describe the physics by MFT and RHA, respectively, 
one fixes the coupling constants differently in both ap- 
proaches. Namely in that way that the minimum of the 
corresponding curve coincides with the phenomenologi- 
cal values fc-O = 1.42 fm"\ {£/pB - to) = -15.75 MeV 
Q. Here, however, our intention is different. Consider- 
ing the Walecka model as underlying rclativistic QFT, 
we see the MFT as ad-hoc description, simply omitting 
the Dirac sea. The RHA, accounting for the Dirac sea, 
is fitted to the phenomenological values of the minimum. 
Thereby the coupling constants are fixed and used for all 
curves (RHA, MFT, NSET). Hence, Fig. Hrather shows, 
how the NSET curve approaches the RHA curve by in- 
cluding higher-order terms in an series expansion in fc/. 

For small Fermi momenta fc/ < 1.1 fm~^ all the curves 
agree nicely but start to deviate for larger values of fc/. 
While the MFT and RHA curves share the same general 
characteristics, namely a local minimum and a steeply 
rising behavior toward larger Fermi momenta, the po- 
sitions of the minima differ by a Fermi momentum of 
= 0.27 fm~^ and an energy = 14 MeV. Including the 
first NSET correction term ~ {'4>'4')^ in. Eq. the cor- 
responding curve valid up to ©(fc'j.^) clearly approaches 
the RHA curve. 
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{£Ipb - m) [MeV] 




FIG. 7: The energy per nucleon in uniform nuclear mat- 
ter with A'' = 2 in 3+1 dimensions determined in full MFT 
and RHA, respectively, is confronted with the NSET results 
valid up to 0{ky) and ©(fc^*). The coupling constants are 
chosen to fit the value and position of the minimum of the 
RHA curve to the phenomenological values k'j = 1.42 fm~^, 
[S / pB—m) — —15.75 MeV denoted by the cross symbol. With 
m = 939 MeV, m„ = 550 MeV, m„ = 783 MeV, we obtain 
= 62.89 and = 79.78. 



C. Spatially nonuniform systems 

Let us extend the analysis of uniform nuclear matter 
in its ground state to include spatially inhomogeneous 
systems. This is necessary to describe finite "nuclei" or 
bound states. 

One should note that it has been shown within 
the Walccka model, that the ground state of infinite 
nuclear matter at a given finite Fermi momentum kf 
does not necessarily imply the fermion condensates to 
be translationally invariant. However, this is not in con- 
flict with the NSET approach, which is based only on 
the assumption of translational invariance of the vacuum 
(kf = 0). By accounting for spatial inhomogeneous con- 
densates in the sector, we rather explicitly allow for 
the formation of such a translationally noninvariant state 
at finite kf. 

While the MFT approach to stationary, spatially uni- 
form matter can be considered as leading-order treatment 
within the NSET approach, the situation is different for 
spatially nonuniform systems. Here, the respective MFT 
equations cannot be inferred directly from the NSET La- 
grangian. In contrast to our NSET Lagrangian, the as- 
sumption of an almost local effective theory, explicitly 
restricting to small momentum transfer, is not yet imple- 
mented in Eq. ([5^ . It is instructive to recall Eq. 




1-^ li' 



FIG. 8: Illustration of the effective fermion current {xpyi^ip). 
a) Valence part {'ilj^^tp) + . b) Residual contribution to the 
fermion current starting with a " — " loop. It is given by 
the infinite sum of topologically distinct diagrams, made up 
of further internal " — " loops and outermost " -1- " loops, 
coupled via a and uj propagators. 

and note that the scalar self-energy can be written as 

Sik) = m + glA'ik){{^^)k + --- . (69) 

Here the ellipses denote further contributions, which are 
of higher powers in the fermion condensates. The vector 
part is determined analogously 

K(fc) = -glDl^^{k){i^^^^)k + ■■■ . (70) 

In the Walccka model the fermion current is conserved 

d^i^-Z^iP) =0 o r(7/;7,i/;)fe = 0. (71) 

Note that ip in Eqs. ((7T|) is the full fermion field. In con- 
trast to uniform nuclear matter, it is not clear a priori 
that an analogous relation for " -I- " states holds for spa- 
tially nonuniform systems also. The expression for the 
fermion current translates into an infinite sum of topolog- 
ically distinct diagrams (cf. Fig.[51). Turning to Eq. ([7T|) . 
the momentum fc" enters the initial " — " loop, Fig. |5)d, 
in the form of fc'^7,y ~ ^. It can be shown that the contri- 
bution of a given " — " loop, with arbitrary 1 and 7^ in- 
sertions, vanishes if one inserts an additional momentum 
via ^. To account for all topologically distinct diagrams, 
one has to sum over all possible insertion points. 
As a result one obtains 

5''(^7.^)+-0. (72) 

For stationary systems, Eq. ([72]) becomes \i'{ipjip)+ = 0. 
This is in particular fulfilled for (7/177/')+ = 0, correspond- 
ing to "positive-energy matter" at rest, which we will con- 
sider in the remainder of this section. 

Transforming into position space and inserting the ex- 
plicit expressions for the effective meson propagators, as- 
suming stationary condensates, we obtain 

Six) = m+gl (^^^ + O ((V2)2)^ (^^7/)+ • • • (73) 

and 

^°(^) = 9l (^T^ + ^ ((V')')) • • • , (74) 

which enter the equation of motion for the " + " fermion 
field 

(i^ + 7oF°(x)-5(x))7/' = 0. (75) 
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Concentrating on Eqs. (|M1) and (US]), MFT dealing 
with stationary, nonuniform systems and the NSET ap- 
proach can be compared. 

The MFT result can be recovered from the NSET ap- 
proach by omitting both the O ((V^)^) terms in Eqs. ([75)1 
and ((74)) , as well as the Dirac sea- induced terms of higher 
loop order, involving increasing powers in the fermion 
condensates and effective propagators, indicated by the 
ellipses. In the framework of the NSET, we do not see 
how such an "approximation" could be justified. 

Let us finally turn to the special case where the meson 
masses dominate the denominators of the meson propa- 
gators. This was assumption in the construction of our 
NSET Lagrangians. Only then do the interactions be- 
come almost local and effective coupling constants can 
be defined. The equation replacing Eq. (|75|) is 

\i0 _ f 1 + ^ + C024^ + • • • ) (VV> 

{ rni \ mi J 

-m+-^ 1 + — + C204^^ + --- W> (76) 

+ ••• I V = 0. 

The ellipses in the last line denote contributions that are 
of higher powers in the fermion condensates, descending 
from effective interactions not present in the original La- 
grangian. Equation (|62p becomes 

|i^ - 470 f 1 + 4 + ^ + • • • ) (V^V.) (77) 
I ml \ ml mt J 

-jn+4 + ^ + ^ + (V5V')U = 0. 
m^ \ mj. mj; / J 

Keeping terms up to ©(V^) in the expansion of the meson 
propagators only, Eqs. ([75)1 and ([77|) coincide. This is a 
consequence of our renormalization procedure. We chose 
a2, cXa and such that up to this order the effective 
meson propagators in the NSET resemble their counter- 
parts in the underlying relativistic field theory. Turning 
to higher-order terms, differences are unveiled. 

The coefficients of the (V^)^ derivatives of the conden- 
sates differ by contributions proportional to the square 
of the coupling constants. Moreover, the effective multi- 
fermion interactions, arising in the construction of the 
NSET, are not accounted for within the MFT. 

Hence, in contrast to the MFT approach to uniform 
matter, in general the MFT approach to nonuniform sys- 
tems cannot be seen as a systematic leading-order treat- 
ment within the framework of the NSET. We believe that 
the NSET approach provides us with a new tool to ac- 
count for the effects of the Dirac sea in a controllable 
manner, avoiding ad hoc approximations. Applications 
of this method to spatially nonuniform systems will re- 
quire numerical studies and have not yet been performed. 



VI. SUMMARY AND CONCLUSIONS 

In this article, we have constructed a no-sea effective 
theory for the Walecka model. We focused on the rel- 
ativistic Hartree approximation. While the Dirac sea 
is integrated out, the positive-energy fermion states are 
treated explicitly. The resulting NSET Lagrangian fea- 
tures an infinite number of effective interaction terms. 
For the case of uniform nuclear matter the NSET could 
be tested quantitatively. Restricting to 1-1-1 dimensions, 
new analytical insights in uniform nuclear matter were 
obtained. 

Finally, we confronted our NSET approach with the 
standard MFT approach with the no-sea approximation 
to the Walecka model, wherein the Dirac sea is simply 
omitted. It turned out that MFT applied to uniform nu- 
clear matter at small Fermi momentum can be reinter- 
preted as leading-order treatment of the NSET. We have 
determined the order in a series expansion in fcy to which 
the MFT results can be trusted; e.g., for the 3-1-1 di- 
mensional Walecka model, the energy per nucleon deter- 
mined by the MFT approach is valid only up to 0{kj^). 
Within the NSET approach, the order of validity can be 
improved systematically. Extending the analysis to spa- 
tially nonuniform systems, the situation turned out to be 
more complicated. Here we emphasized the significant 
differences between the two approaches. 
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Appendix A: Exact localized solutions of no-sea 
effective theory with four-fermion-interactions 

In this section we use an approach, introduced by Lee 
et al. [l3 to construct exact localized "classical" solu- 
tions of 1+1 dimensional field theories with four-fermion- 
interactions. To apply their considerations to Eq. (PH]). 
we have to generalize their ideas to the simultaneous ap- 
pearance of scalar and vector interactions featuring two 
different coupling constants Ga ^ G^^, 

£ = (i^ - m) V + ^ {i^i^f 

+ ^ (^i7,V) (Vii7» • (Al) 

Here "0 corresponds to a " -|- " fermion spinor with N 
fiavors. The required generalization is straight-forward. 
Following Lee et al, fermion bound states can be con- 
structed. To retain analyticity, they assume that all 
fermions occupy one and the same energy level (energy 
eigenvalue E). Obviously, this limits their approach to 
fermion bound states with n < N fermions only. It is 
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convenient to turn to the notation of [l^l- This imphes 
the following choice of the 7 matrices 

7o = CTg , 7i = -i(Ti , 75 = 7o7i = 02 , (A2) 

and the ansatz (fc denotes the flavor index) 

^7ifc\ R(x) ( cos ff^ 



Vk 



/cos e\ 



(A3) 



where Ufc, Vk have to be real valued functions ([T3|, below 
Eq. (7c)). 9{x) is determined by 



6{x) = arctan (a tanh/3a;) 



(A4) 



a and /? are energy- and mass-dependent constants de- 
fined by 



rn — E 
m + E 



and /3 = \/r 



E^. 



(A5) 



By using Ref. [l^, Eq. (8a), together with the La- 
grangian (jAip . we infer 



2 (£; - m cos 26*) 



G2 - Gl cos2 26* 
Demanding the following normalization 



(A6) 



n — dx i/j'^iIj 



de 



G2, - Gl cos2 261 



, (A7) 



where n > 0, Eq. (jA7p can be fulfilled consistently only 
for G?. > G^,. In this case we obtain 



2 arctana = arctan 
where we defined 



^tanh (|g^G 



G=VGl-Gl. 
Equation ()A8|) can be also written as 
a G 



1 - a2 2 G^ 



■ tanh 



(|g„g) 



(AS) 



(A9) 



(AlO) 



The "bound-state mass" of n fermions corresponds to the 
expectation value of the classical Hamiltonian. It is given 
by (cf. Ref. [3, Eq. (7c)) 



arctan a 

Af„ = m I dx tpip = m / d9 



-2 cos 26* 



G2 - Gl cos2 261 



2 m 
G(j G 
2m 



arctanh 
arcsinh 



— ^ sin (2 arctan a) 
G 



Ga 

Gu. 



sinh 



(All) 
(A12) 



In the last step of Eq. (|A12|) we inserted Eq. Let 
us investigate the stability of the n-fermion localized so- 
lution. To be stable, it has to be energetically favored as 
compared to the mass of n free fermions, as well as the 
mass of n localized fermions not interacting with each 
other. The first condition is true, as 



Arsinh(w sinh z) < 



(A13) 



for w > 1 , z > directly implies M„ < nm. To fulfill 
the latter condition, we have to ensure that M„ < nMi. 
This is also true as 

Arsinh(w sinh(7iz)) < 71 Arsinh(w sinh z) (A14) 

for w > 1 , z > 0. In result, the n-fermion localized 
solution with mass Af„ is energetically favored. Finally, 
wc determine the energy eigenvalue E. It is given by 



E = m 



I -a' 
l + a2 



(A15) 



Recasting Eq. (|A11|) into 
a 1 



2 — 2 sin (2 arctana) 



2Ga V 2m 
and using Eq. (|A10|) . one obtains 



^ tanh(^^Af„) (A16) 



cosh[%G^^Gl-Gl 



l+(f^) sinh^i^f G^yGr~G2 



(A17) 



Therewith, the explicit expression for the n-fermion-state 
spinor becomes 



ip{x) 



R{x) 



\/l + a2 tanh^ (3x V" tanh (ixj ' 



1 



(A18) 



where we insert a and fi defined in Eq. (jASP with E 
from Eq. (|A17[) . Note again that to be a solution, ip{x) 
has to be real valued. Adopting the explicit values for 
the coupling constants, we obtain a solution of the La- 
grangian (j40p . we arc interested in here. For the special 
case Ga = G^ we refer to Ref. p^ . 



Appendix B: Results for coefficients defined in 
Sec.lrVl 



Here we collect the coefficients associated with the n- 
fermion bound state with mass Af„ in Sec. IIVI In order 
to simplify the notation, we define 



A 



Agl 



B = 



Bgl 



c^A-B. 



(Bl) 
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CocfRcicnts for the scalar potential S{x), Eq. (15^ . 
n'^Ac n'^Ac^ 



S22 



S44 



S52 



S54 



4 



S42 



96 



55) , 



64 
n^Ac 
72 



(i + 35) , 



96 



TT 
TT 



4AA - iAB ~ BB 



(35 - 5i) 



^''"'^ ^5Ai - 3A5 - 255^ 



S62 
S66 



23040 



3072 

n'^Ac^ 
1024 



135^^ - 210yl5 + 915' 



9^^ - 14^5 - 275' 



^2 + 10A5 + 552 



(B2) 



CocfRcients for the potential Vb(x), Eq. (|57p . 



P22 

Pi2 
PiA 
P52 

P54 = 

P62 = 
P64 = 
P66 = 



n?Bc 
n^B^c^ 



48 

64 

72~ 

n^Bc 
l8~ 



(3^ + 5) , 



AA + SAB - ABB 



NA^ Gm'^c / , 7 7 
+ — ^ 2AA + 3AB - 5BB 



n^B^c^ 

1440 
n^Bc^ 
' 3072 



15^-^ - 42A5 - 55' 



n^Bc^ 
1024 



5^2 + 10A5 + 5' 



(B3) 
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